Detailed uncertainty analysis of the ZEM-3 

measurement system 
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Objectives 

Develop an uncertainty 
analysis for a common 
resistivity and Seebeck 
coefficient measurement 
configuration. 

• Introduce a software package 
which includes the 
uncertainty analysis 
calculations. 

•Establish measurement best 
practices to minimize 
measurement uncertainty. 
Demonstrate typical high 
temperature uncertainty on a 
Si/Ge sample. 



Potentiometric 
Configuration (4-probe) 

Power Factor Uncertainty 
+ 7 % / - 25 % 
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Resistivity Uncertainty 


Sample Uniformity 


# 

Source 

Typical Values 

1 

Thermocouple tip radius 

0.25 mm 

2 

Thermocouple separation length 

±0.1 mm 

3 

Sample uniformity 

±0.1 mm/cm 

4 

Caliper uncertainty 

±0.01 mm 

5 

Statistical variation 

Calculated 

6 

Wire discrepancy 

Calculated 

7 

DAQ voltage uncertainty 

50 ppm +1.2 pV 

8 

DAQ current uncertainty 

0.2% +0.3 mA 



Typical Thermocouple Beads 



Source: ULVAC Source: Cleveland Electric 
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Seebeck Uncertainty Sources 


Cold-finger Effect 


# 

Source 

Typical Values 

1 

Cold-finger effect 

10,000 W/(m 2 K) 

2 

Wire Seebeck variation 

±5% 

3 

Absolute temperature 

±2°C 

4 

Statistical variation 

Calculated 

5 

Wire discrepancy 

Calculated 

6 

DAQ voltage uncert. 

50 ppm + 1.2 pV 

7 

DAQ temperature uncert. 

50 ppm + 1.2 pV 


Wire Seebeck 




Upper Probe 


|^3 Lower Probe 


Close-up View 
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Resistivity Calculation 

_ 2 Zi 2 Vi ~ N Y. z iYi wD 
P ~ (S Zf ) 2 — NYjZ t 2 L 

y- probe to probe voltage 
z- test electrical current 
w/D- sample width/depth 
L- probe separation length 
N- sample size 


Seebeck Calculation 


-NY.Xiyt o 

x- probe AT 

y- probe to probe voltage 
S W ire" w ' re Seebeck coefficient 
T- sample temperature 
N- sample size 


Error Propagation 



Uncertainty can be calculated from a 
Taylor Series expansion around the 
nominal measurement value 
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FEA Model Parameters 

•Autodesk Simulation Multiphysics 

• Non-linear iterative thermal solver 
•Thermal domain: 

•Sample 4x4x18mm rectangular prism 
•Thermal conductivity 4 W/(m 2 K) 
•Probes 00.5mm x 150mm 

•Thermal conductivity 30 W/(m 2 K) 
•Boundary conditions: 

•Sample ends fixed temperatures 
•Probe ends fixed temperatures 
•Remaining faces radiation coupled 
•e=0.2, 0.5, 0.7, 0.9 (40% change) 

• Parameters of study: 

•Furnace temp=200,600,1000°C 
•Differential temp=0.1 to 14°C 
•Thermal conductance= 100,000, 33,000, 
10,000 W/(m 2 K) (600% change) 


Grid Independence Study 



•Two meshes were generated 
from primarily brick elements 
•Course mesh (shown above) 
•41,000 elements 
•Fine mesh 

•55,000 elements 
•Mesh agreement <0.2% change 
in results 
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Temperature Contour 


Sample 


wm 1006.0 C 
H- 1004.8 
1003.6 
1002.4 
1001.2 
1000.0 


Upper Probe 


Iti T2 


Lower Probe 


• T1/T2 represent "actual" and 
"measured" temperatures 
•Model fits experiment well at 
high temperature 


FEA Uncertainty Results 

•Uncertainty is defined as the difference 
between the desired and measured 
temperature difference. 

•Uncertainty increases with furnace 
temperature and delta temperature. 
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Delta Temperature 
Uncertainty (Celsius) 



ise 


Testing Samples & Profile 

Si 80 G e 2 o samples prepared by milling and 
spark plasma sintering elemental 
powders 

•2at% P doped 

•0 1" pucks machined to 4x4x18mm 
Samples measured from 25 to 950°C 
Equilibrium definition: 

•Furnace <5% change in 120 seconds 
•Isothermal <0.1°C in 120 seconds 
Resistivity measurement: 

• -50 to +50 mA increment 5mA 
Seebeck measurement: 

•+l°C/min up to 10°C 


Example Measurement Data 




Error in Resistivity (%) 


ise Study 


Resistivity Uncertainty 



0.5% Error 


Resistivity uncertainty is fairly 
temperature independent, due to 
geometric nature. 

Thermocouple tip radius dominates 
uncertainty. 


Seebeck Uncertainty 



- - ♦ 0.5% Error 


Seebeck uncertainty is highly 
temperature dependent. 
Cold-finger effect dominates at all 
temperatures, and is asymmetric. 
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Resistivity (Ohm-cm) 


Oo^irodlMeeSoirD S®Mir©s§ IPir®p®$i)ftii®irD (FIM IRtodteltog Case Study 


Overall Results 


Power Factor Results 



CO 


Resistivity uncertainty is ±7.0% at all 
temperatures. 

Absolute Seebeck uncertainty ranges 
from ±1.0% at room temperature to 
+1.0%/-13.1% at high temperature. 



* Power factor uncertainty ranges from 
±7.5% at room temperature to 
+7.3%/-25.0% percent at high 
temperature. 

These values all assume the 
conservative parameter values listed. 
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Conclusion 

LabVIEW VI is available to operate the 
ZEM-3 and calculate the uncertainty. 
Resistivity uncertainty is primarily 
geometric, and can be reduced with 
careful preparation. 

Seebeck uncertainty is primarily due to 
the cold-finger effect, and can be reduced 
with good thermal contact. 

Power factor uncertainty is ±7.5% at room 
temperature and +7.3%/-25.0% percent at 
high temperature. 
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Error Propagation 


Resistivity and Seebeck 


y . 



^ 'Total 




+ ... 


E Zi E Vi ~ N T. z t yi wD 
(Zz.y-NZzt 2 L 


S = 


Zx t 'Zy i -N'Zx t y i o 
( ZXiy-NXxy + *™reU) 


x- probe temperature difference 
y- probe voltage 
z- electrical current 
N- sample size 


Statistical Uncertainty 

Ustat ~ ^v, 95% 

NI(y; -y c 0;)) 2 

v(WSz, 2 -(Szi) 2 ) 
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Dimensionless Temperature (T/Tmax) 


FEA Modeling €a§@ §ta(% 


FEA Verification 



Temperature 

Model 

Experiment 

Celsius 

ATprobe/AT 

ATprobeM"T 

200 

0.33 

0.18 

600 

0.17 

0.10 

1000 

0.07 

0.06 


Temperature Contour 

T1 desired temperature 
T2 measured temperature 
• accounts for thermal 
contact conductance 



1 1006.0 c 

1004.8 

1003.6 

1002.4 

1001.2 

1000.0 

Upper Probe 


Sample 


Lower Probe 
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Furnace Test Temperautre (Celsius) 



FEA Uncertainty Results 

•Uncertainty is defined as the difference 
between the desired and measured 
temperature difference. 

•Uncertainty increases with furnace 
temperature and delta temperature. 


0.21 


0.55 


Probe Delta Temperature (Celsius) 


1000 


Influence of Thermal Contact 
Conductance 

•Thermal contact conductance 
plays a significant role in the 
Cold-finger effect and displays a 
power law dependence. 



• 1 1 1 1 

0 20000 40000 60000 80000 100000 

Thermal Contact Conductance (W/(m 2 K)) 

Furnace Temperature 1000°C 
Delta Temperature 14°C 
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